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^ ' Abstract 

This paper promotes the differential method as a new fruitful strategy 
for estimating a ground-state energy of a many-body system. The case 
. of an arbitrary number of attractive Coulombian particles is specifically 

' studied and we make some favorable comparison of the differential method 

, to the existing approaches that rely on variational principles. A bird's-eye 

■ view of the treatment of more general interactions is also given. 

PACS : 24.10. Cn, 12.39. Jh, 03.65.Db, 05.30.Jp 

^ ■ 1 Introduction 

There is little need to stress the great importance of the ground-state in many 
domains of quantum physics. Nevertheless, computing the lowest energy of 
' most systems cannot be done analytically and approximations are required. 

M I Some techniques, being very general, have been well-known for several decades. 

For instance the variational methods (Rayleigh-Ritz) or the perturbative series 
(Rayleigh-Schrodinger) have their roots in the pre-quantum era; much later, 
numerical algorithms (numerical diagonalizations as well as Monte-Carlo com- 
putations), supported by the increasing power of computers, have been able to 
provide a tremendous precision on the ground-state of a large variety of very 
complex systems. However, it is a much more difficult task to rigorously es- 
timate the discrepancies between the exact ground-state energy Eq and the 
approximated one. In particular, though variational methods naturally pro- 
vide upper bounds on Eq, obtaining lower estimates requires more sophisti- 
cated techniques (for instance the Temple-like methods [191 § XIII. 2]), some of 
them being very system-dependent {e.g. the moment method proposed in |1(J | 
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for rational-fraction potentials or the Riccati-Pade method proposed in f7| for 
one-dimensional Schrodinger equations). 

The optimized variational methods. For a many-body system governed 
by pairwise interactions, an interesting strategy is to approximate Eq from below 
in terms of the ground-states of the two-body subsystems . Such an approach 
has been successfully applied to Coulombian (bosonic and fermionic) systems of 
charged particles or self-gravitating bosons ^ . Clever refinements have 
been proposed that provide some very accurate lower bounds of Eq for the three- 
body 13 and the four-body systems Though not easily generalizable to an 
arbitrary number of particles, these last optimized variational methods can be 
applied to interactions that are not necessarily Coulombian and may be relevant 
for quarks models [2] where some inequalities between baryon and meson masses 
represent theoretical, numerical and experimental substantial information |17[ 
1201 and references therein] . In practice, the optimized variational methods allow 
to efficiently treat some models that have simple scaling properties, for instance 
when the two-body interaction can be described by a purely radial potential of 
the form v{r) oc sign{f3)r^ . The main reason relies in the fact that, except the 
Coulombian {(3 = —1) and the harmonic (/3 = 2) interactions, the exact form 
of the ground-state energy of the two-body problem is not known. Yet, one can 
still take advantage of the power-law behavior of v to obtain worthwhile lower 
bounds for the ratio between the A^-body and the 2-body ground-state energies. 

The differential method. Besides the variational and perturbative tech- 
niques that are mentioned above, there exists a third very general method for 
approximating the ground-state energy of a quantum system, namely the differ- 
ential method (see El and references therein for a historical track] ) whose 
starting point is recalled in sectional for the sake of completeness. As for the 
variational methods, the differential method call on a family of trial functions 
that supposedly mimic the ground-state and that allow for the construction of 
a function (the average of the Hamiltonian in the former case, the so-called 
local energy in the latter case) whose absolute extrema within the chosen trial 
family provide bounds on the exact ground-state energy Eq. One of the main 
advantages of the differential method over the other ones is that no integral is 
required and, then it allows to work, even analytically, with rather complicated 
trial functions by encapsulating some rich structure of the potential. It is also 
worth mentioning that the same test function leads to both upper and lower 
bounds on Eq and then the estimates comes with a rigorous window. Though 
applicable to many models (to systems involving a magnetic field, to discrete 
systems, to non-Schrodinger equations, etc.) the major inconvenient of the dif- 
ferential method is that it requires, as a crucial hypothesis, the exact eigenfunc- 
tion to remain non- negative in the configuration space (it must have no node). 
Therefore it excludes any ground-state whose spatial wave-function is antisym- 
metric under some permutations of its arguments. As far as fermionic systems 
are involved, the differential method will concern only those whose ground-state 
eigenfunction remains symmetric under permutations of the spatial positions of 
the identical particles. 

The aim of this paper is to apply the differential method specifically to a 
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system made of A'' non-relativistic particles of masses rrii in a D-dimensional 
space whose Hamiltonian has the form 




(1) 



When the N particles located at {ri}i=o,...,A'-i interact only through pairwise 
potentials Vij = Vji , V is given by 



where = rj — r^. The spin-dependent interactions, if any, are assumed to 
be included somehow in the scalar potential V and H will be supposed to act 
on spatial wave-functions only; in other words the possible spin configuration 
have been factorized out in one way or another. When iV = 3 and = 4 and 
for power-law potentials, this is the same kind of systems to which the opti- 
mized variational method applies also. We will consider the Coulombian case 
in section and systematically compare the estimates given by the variational 
methods and the differential method. In section 01 general interactions are 
considered (not necessarily power-law Vij 's) . This is of course relevant for esti- 
mating the ground-state energy of a system where the spin-independent strong 
interactions are dominant; for heavy enough quarks for instance, it is known 
[TT] that the non-relativistic form ^ may be pertinent^. At atomic scales, the 
method could be applied to clouds made of neutral atoms where short-range 
interactions govern the dynamical properties. 

2 The differential method 

2.1 The general strategy 

The necessary but sufficient condition for the differential method to work is 
the following: the Hamiltonian H has one bound state |$o): associated with 
energy Eq, such that ^o{q) = ('j'I'I'o) remains non-negative in an appropriate 
g-representation, say of spatial positions. For a A^-body system governed by 
the Hamiltonian J^l, the dynamics in the center-of-mass frame corresponds to 
a reduced Hamiltonian H whose ground-state^ $o has precisely this positivity 

^Possible relativistic corrections may be included (for instance by considering the spinless 
Salpeter equation) since the differential method does not require a quadratic kinetic energy. 

■^We will only consider the cases where at least one bound state exists. Physically, this 
can be achieved with a confining external potential (a "trap" is currently used in experiments 
involving cold atoms). Formally, this can be obtained in the limit of one mass, say mo, being 
much larger than the others. The external potential appears to be the VQi^s, created by such 
an infinitely massive motionless device. It will trap the remaining N — 1 particles in some 
bound states if the VQi's increase sufficiently rapidly with the roi's. 
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i,j=0 
i<3 
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property in the whole configuration space Qat of the (A^ — 1)d relative coordi- 
nates gjv '= (i"! — To, . . . ,rjv-i — To). This is the Krein-Rutman theorem (see 
[l9. §XIII.12]). For each state \lp), the hermiticity of H implies the identity 
($ol(^ ^ EQ)\^p) = 0. If we choose \(p) such that its representation ip{qN) is a 
smooth normalizable real wave-function, we obtain 

/ <^l{qN){H-E^)<p{qN)dqN^Q. (3) 

Jqn 

Taking into account the positivity of $o on Qn, there necessarily exists some 
qN such that {H — Eo)(p{qj\i) ^ and some other configurations for which 
{H — Eo)^p{qN) ^ 0. Choosing </? > on Qn, we get both an upper and a lower 
bound on Eq: 

inf {El:liqN)) ^ Eo ^ sup {El:liqN)) , (4) 

QjV Q„ 

where the local energy is defined by 

. dof HipjqN) , , 



In other words, the differential method provides an estimate 



E, 



(d.m.) dcf 
~ 



.Qn 2" 



that comes with a rigorous windows ±A_Bq'^ '" '' where 



AE, 



(d.m.) dcf 
~ 



SUp{El:liqN))-M{El:liqN)) 
.Qn 2" 



(6) 



(7) 



Unlike for the variational method, the determination of the absolute extrema 
of the local energy does not require the computation of any integral. Even the 
norm of the test function ip is not required provided it remains finite. The 
two inequalities Q become equalities (the local energy becomes a flat function) 
when 93 = $0 and therefore we will try to construct a test function that mimics 
$0 at best. We will choose ip that respects the a priori known properties of 
$0: its positivity, its boundary conditions and its symmetries if there are any. 
Since for each test function the error on Eq is controlled by inequalities Q), the 
strategy for obtaining decent approximations is clear: First, we must choose 
or construct (p to eliminate all the singularities of the local energy in order to 
work with a bounded function. For instance, when the Hamiltonian has the 
form + V with V being unbounded at some finite or infinite distances, the 
first kinetic term of the local energy i^l^j = —A(p/(p + V must compensate 
the singular behavior of V for the corresponding configurations (we will work 
systematically with units such that fi = 1). Once a bounded local energy, say 
eI^^\ is obtained, we can proceed to a second step: perturb the test function, 
fQ—^^ — ^o+ '5(y3, in the neighborhood of the absolute minimum (resp. maxi- 
mum) of eI^°^ in order to increase minE'l^^."' (resp. decrease maxE'l^^."'). Up to 
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the end of this article, we will focus on the first step: we will show how obtaining 
a bounded local energy furnishes some sufficiently constrained guidelines for ob- 
taining reasonable bounds on Eq'^. We will keep for future work the systematic 
local improvements of the absolute extrema of the local energy. In [151 § 6], I 
have shown on a simple example how this can be done. 

2.2 Illustration in the two-body case 

Before coping with much complex systems, let us first consider the case of the 
two-body problem that can be reducible to a one single non-relativistic particle 
of unit mass in an external potential V. The local energy is 

4:1 --^ + ^ = - i(V5)^ + V (8) 

where S '= ln((p) is a well-defined function when ip > 0. 

For the D Coulombian potential V{r) — n/r (d > 1), the singularity at r = 
controls the local behavior of the test function if one wants a bounded local 
energy. If lim^^o S{r) is finite, by possibly subtracting an irrelevant constant 
term, we can suppose that this limit vanishes. Therefore, without too much loss 
of generality, we assume that 5, can be asymptotically expanded near r = 
on a family of power functions (for D = 1, the logarithmic functions should be 
considered) whose dominant term can be written like S{r) ^^r^o s^r'^'^^ / {a+\) 
with fT 7^ — 1 and so 7^ 0. Balancing the dominant terms in the local energy, it 
is therefore straightforward, to check that the only choice for the parameters so 
and a to get rid of the Coulombian singularity is too take a = Q and so = 
2k/ {u — 1). It happens that for S exactly equal to 2Kr/{n — 1), we obtain a 
global constant local energy, namely — 2k^/(d— 1)^. Therefore we have obtained 
the exact wave-function of the ground-state provided we eventually check that 
the wave- function is square-integrable which is true for k < 0. 

For the harmonic oscillator V{r) = V is unbounded as r increases 

(r = +00 is a singular point for V). If we tentatively look for an S whose 
asymptotic expansion at r — s- -\-qo has a leading term of the form sor'^+^/((T + l), 
we necessarily get a = 1 and so = — w (the cases where cr ^ — 1 are ruled out 
by the square-integrable property). The local energy Da;/2 is actually constant 
for all r's and indeed represents the exact ground-state energy. More generally, 
with the help of standard linear algebra arguments, for V being any definite 
positive quadratic form, we can always find a quadratic form S for which the 
local energy is globally constant. 

^One can understand it from the extreme sensitivity of the local energy to any local per- 
turbation of the test function: while, in the variational methods, the quantity {ip\H\if) / {ip\^p) 
is quite robust to local perturbations because it represents precisely an average on the con- 
figuration space, the local energy may become unbounded quite easily by canceling locally f 
faster than Hif. 
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2.3 Formulation for the many-body problem 

When the potential V has the form a natural choice of trial function is to 
take (for variational techniques in a few nuclear body context, such a choice has 
been used by for instance) 

l<j 

where each of the N{N — l)/2 functions 4)ij{v) — (f>ji{—r) depends on D coordi- 
nates^. 

It is straightforward to check that this choice describes a state with a fixed 
iter-of-mass: indeed v 
local energy is given by 



center-of-mass: indeed we have {J2iLo^ Pi)\v) = 0- Hence, H^p = Hip and the 



N-l 



E 



(10) 

1 V(/)y(ry) V0jfc(r,fc) 



where m.y stands for the reduced masses mimj/lmi+mj). The last sum involves 

all the N{N~ 1)(A^— 2) /2 angles {j, i, k) between and r^^ that can be formed 
with all the triangles made of three particles having three distinct labels (z ^ j, 
i ¥^ k, j k). Let us now take to be a positive solution of the two-body 
spectral equation 

-^^^'/'y + = ty'/'^i • (11) 

The local energy becomes 

eI%n) = Yl - E — ^S^,{r,,) ■ VS.,k{r,k) (12) 

- ■ „ - - 

where 5*^^ ln(0ij). The trial wave- function Q of the global system must 
be kept square-integrable but it is not necessary for all two-body subsystems 
to have a bound state when isolated^. For instance, for two electric charges 

*The present paper wants mainly to stress the simplicity of the differential method. It 
does not seek for a real performance at the moment and we will not try to improve the choice 
of coordinates. Working with Jacobi coordinates, for instance, or constructing optimized 
coordinates as done in ^ may lead to better results. Anyway, we will see that the numerical 
results of section |3 are satisfactory enough for validating the approach by the differential 
method. 

^But for some pairing, 1111 must have a normalizable solution. The cases of Borromean 
states where no two-body binding is possible l21l cannot be described by the form if we 
keep llli . 
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having the same sign a positive but non-normahzable solution of Hll|l can be 
found. When Vij admits at least one bound state (see also footnote thanks 
to the Krein-Rutman theorem, we are certain to get a positive (jjij when taking 
the ground-state of the two-body system and its corresponding energy. At 
finite distances, if the possible singularities of Vij are not too strong, we expect 
that 4'ij and then Sij to be smooth enough for Eioc to remain bounded. At 
infinite distances, -Eioc is expected to become infinite if Vij does not tend to a 
constant sufficiently quickly. To see that, one can take purely radial potentials. 
i.e Vij{r) = Vij(r) where r =^ ||r||, and consider the asymptotic behavior of Sij 
given by the semiclassical (jwkb) theory (see for instance J^). Its derivative is 
given by S'ij(r) ~r-^cx3 — \/2my [?;ij(r) — e^] and is not bounded if Vij is not (at 
infinite distances). Strictly speaking, it is only for short-distant potentials that 
we can hopefully obtain rigorous non-trivial inequalities (@J while keeping the 
choice © with f(T^\ . However, as will be discussed in section Q), the ground- 
state energy may be generally not be very sensitive to the potential at large 
distances (far away where <I>o is localized) and this physical assumption may be 
implemented by introducing a cut-off length from the beginning. 
For purely radial potentials H12|l simplifies in 

^ in,)Sl,ink) cosilXk) . (13) 

_ — ... mi 
j,i,k 

Yet, for a multidimensional, non-separable, Schrodinger equation like 
a JWKB-like asymptotic expression is generally not available |14[ Introduction]. 
Nevertheless, the differential method is less demanding than the semiclassical 
approximations: we will try to keep the local energy, like the one given by (|12|l . 
bounded at infinity but we will not necessarily require it to tend to the same 
limit in all directions. 



Af-l 



3 The Coulombian problem 

The purely Coulombian problem in D > 1 dimensions corresponds to the situa- 
tion where all Vij 's are radial potentials and have the form 

v.,{r) = '-f (14) 

for N{N — l)/2 coupling constants that may be or may be not constructed 
from individual quantities like charges. Provided a A^-body ground-state exists, 
we can solve exactly Hll|l making use of A0(r) — 4)"{r) -I- (d — l)0'(r)/r. We 
obtain a bounded local energy given by 

N-l o 2 , 

- E - 7^ E "^^^^^^cosijXk) . (15) 
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For obtaining upper and lower bounds on Eg, one has just to calculate the ab- 
solute extrema of such a function. It can be done by standard optimization 
routines up to quite large A'' and even analytically in some cases (see below). 
The recipe is therefore simple and systematic: as far as only Coulombian inter- 
actions are involved, we can work with generic masses and coupling constants 
for which @ is normalizable. The remaining of this section will concern the 
quality of these bounds and then we will accord our attention to cases that have 
been treated by other methods, mainly those treated in the references cited in 
the second paragraph of the introduction. More specifically, in order to leave 
aside the problem of the existence of a ground-state we will consider the case of 
attractive interactions only^ (all e^'s being negative). 

3.1 Arbitrary number of identical attractive particles 

In this section we consider one species of particles only: for all i and j we 
denote rrii — m and = —.9^. The local energy H15|) becomes 

E^In) - - Qa^(^ - 1) + F^iq^)^ (16) 

where eo =^ mg^. The function 

FN{qN) J2 cos(jXfc) (17) 

is invariant under translations and rotations but also under dilations of the 
particle configuration. For = 3, the appendix|31proofs that supg^ F3 = 3/2 is 
reached when the three particles make an equilateral triangle and infgj F3 = 1 
is obtained when they are aligned. From this last result we are able to provide 
the lower bounds for for any N: by a decomposition of Fn into a sum 
on N{N - 1){N - 2)/6 triangle contributions, 

-Fjv('?jv) = COs(i3,Zi,Z2) -l-cos(ji,Z2,i3) -l-cos(j2,i3,ii) , (18) 

— ^ ^ ^ 

{11,12,13} _ r / ^ 

i!iii<i2<i3<N — ^3(,r.jj,ri2,ri3j 

all the i^3's in the sum reach their minimum simultaneously when all the particles 
are aligned, and for this configuration we have infg„ F/v — N{N — 1){N — 2)/6. 

^For an immediate application in the case of charged electric particles see IISI where the 
Helium atom is discussed. While the differential method provides an analytical non-trivial 
upper bound for the ground-state of the Helium-like atoms, Eq ^ —{Z — l/2)'^ in the simplest 
model (d = 3, non-relativistic, spinless and with an infinitely massive nucleus of charge Z in 
atomic units), in the case of more delicate systems like the positronium ion (e+,e~,e~), a 
systematic improvement is clearly required but is beyond the scope of this paper as explained 
at the end of subsection l2.1l Indeed, for (e+,e~,e~) if we content ourselves with eliminating 
the singularities, we obtain, for D = 3, —9ma^/8 ^ Eq ^ (o being the fine structure 
constant). The upper bound is trivial while the lower bound is even worse compared to 
—3ma^/4 obtained in eq. (6.4)] or to —Sma^/A obtained by a simple crude argument |21 
eq. (6.6)] (the exact result is Eq < — ma^/4). 
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From we deduce that for each N 

Eo^-^^^^N{N~l){N+l) . (19) 

For D = 3, the same exponential test-functions lead to the better variational 
estimate eq. (17)] 

Eoii-ef,— N{N-l)^ (20) 

(25/512 ~ 0.0488 > 1/24 ~ 0.0417). This was expected from the general 
identity valid for any normalized function (p, 

^*iq)H^{q)dq^ [ |(^(g) psj^^l (g) < sup (i^L^' (g)) . (21) 

The differential method always gives worse upper bounds than the variational 
method with the same test-functions but, in the last case, one still has to be able 
to compute the integrals and one cannot generally estimate how far from the 
exact value the average Hamiltonian is. For a different choice of test functions, 
a better variational upper-bound has been obtained eq. (16)], 

Eq < -.0542iV(iV - 1)2 . (22) 

As far as lower estimates are concerned, bounding from above will allow 
us to improve the existing results, namely (for D = 3) 

Eo^^-^N^iN^l) (23) 

obtained in P eq. (12)]. 

First, when N is not too large for the numerical computation to remain 
tractable, the direct calculation of supg^^ Fn shows (see figure that it gives 
better lower estimates than (|23ll . For very large N, we can nevertheless ben- 
efit from the maximum of Fm for smaller M. Indeed, for M ^ N we can 
decompose Fn into contributions of M-clusters as follows: 

, , (M- 3)! (TV- A/)! , , 
FNiqN)= ^ ^ ^^^^y^ '-FM{qM) (24) 

_A/— subclustcrs 

where the sum is taken on all the Af-subclusters, labeled by the coordinates qi\i, 
that can be formed with the given configuration q^. This sum involves exactly 
N\/M\/{N - M)\ terms and we have 

N{N ~1){N -2) , , 



This leads to define 



dot supg^^ Fm , , 

^ M(M-l)(Af-2) 
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5 10 ^5 N 20 25 30 

Figure 1: Estimations of ground-state energy for N identical attractive 
Coulombian particles (rescaled by a factor N^). The differential method pro- 
vides E'q'^ ™'^ defined by © (black disks) with exact error bars whose half- 
width are given by {Tj). Each solid line starting at M corresponds to the lower 
bounds H27|l . For M = 3 and M — 4, there is only one solid curve given by (|23|) . 
The dashed line corresponds to the upper bound H20|l given by ^| eq. (17)]. 
The dot-dashed line corresponds to the upper bound (|22(l given by eq. (16)] 
that tends, when N — > oo, to the thin dotted horizontal line located at —.0542. 

and from H16|l we find 

^ " (^^ " + ""''^^ " ^0 ■ ^^^^ 

Since, from H25|l . ajv/ is decreasing when M increases, the larger M the better 
the lower estimate of Eq. 

For M = 3 we have already seen that — 1/4; for D = 3, H27|) repro- 
duces exactly (1^ . No better estimate is obtained when considering M — 4. 
Indeed, the configuration of particles that maximizes F4 corresponds to the reg- 
ular tetrahedron because its faces, that are equilateral triangles, maximize the 
contributions of all the 3-subclusters simultaneously. We obtain supg^ F4 — 6 
and hence = a^. 
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For M = 5,6, 7,8 and D = 3, the configurations that maximize Fm can be 
seen in figure |21 Crossed numerics and analytical studies lead to very plausi- 
ble conjectures on the geometrical description of the configuration for M = 5 
and M = 8 for which explicit analytical value of Um can be proposed ^Hl ■ The 
lower bound for (|27|) . for ^ M, is strictly improved when increasing M from 5 
and in particular is better than H23|) . 



N=3 



N=4 



N=5 



A 4^ 



N=6 



N=7 





Figure 2: The configurations that maximize F/v defined by H17|l for identical 
Coulombian attractive particles. 

However the sequence of improvements obtained this way seems to saturate 
up to aoo = 2/9: 

(23 = a4 = ^ ^ as ~ .2432 ^ ae ~ .2414 ^ aj ~ .2382 ^ ag ~ .2366 ^ • • • 

2 

• • • ^ ago ~ .2266 ■ ■ ■ ^ = - ~ .2222 . (28) 

9 

From the optimized configuration obtained with N about several tens (see fig- 
ure |3^), we can hopefully guess that the limit N ^ oo leads to a continuous 
and uniform distribution of the particles on the same sphere. The continuous 
limit of supg^ F/v varies as with N. If, on the unit sphere S, the N par- 
ticles get distributed uniformly with density a = N/Att, the continuous limit 
of supg^ -FV is given by [dSi is an infinitesimal portion of the sphere near the 
point Pi, i = 0, 1, 2; see figurelHb) 
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Figure 3: For large N, is maximized for a configuration where the N points 
seems, numerically, to distribute uniformly on a sphere (with the well-know 
caveat concerning the ambiguity of such a notion Such an optimal config- 

uration for TV = 30 is shown in a) . The continuous limit of -Fjv /N^ when N —> oo 
can be computed with the help of figure b). 



sup Fat - [ adSo^ [ adSi [ a dS2 cos(Pi7^P2) ; 

X /'27r ^ ^ 

-a^ / dai / da2 / d?/'2 sin(ai) sin(a2) cos(Fi, Pq, -f2)- 
2 Jo Jo Jo 

When cos(Pi, Po, ^2) — np°p^| • |p°p^| is expressed as a function of ai, a2 
and tp2, we straightforwardly get supg^ Fn/N^ ^n^oo 2/9 — aoo- For infinite 
TV, this last result supplement the upper bound given by H22|) and we have 

1 E 
N^oo; -— ~-.0556 < ^ < -.0542 . (29) 

3.2 Three particles with one different from the two others 

Without loss of generality, we can choose units where ei2 — 613 = 623 = — 1, 
mi = 7712 — 1, mg = m. For D = 3 the local energy p5|) simplifies into 

4(m + 1) 2(m + 1) \ m + 1 / 
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The general study of appendix [S] applied for (oi, 02, 03) = (1, 1, -p^^) allows to 
get the following analytic bounds on Eq: 

Os^TO=^l; rrf !^Eo^ r ^ ; (31a) 

' 4 8 (m + l)2 ^ " 4 (to + 1)2 ' ^ 

1 m m(m + 2) 1 mhn + 2) , , 

1 ^ m < 3; - V / , ^ ^0 -T - , ,,2 ; (31b) 
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3 < to; 7^ TTT^ ^Eo !^ H • (31c) 

' 4 (to + 1)2 ^ 4 (m + l)2 ^ ' 

For TO ^ 3, the lower bound on Eq corresponds to a configuration where the 
particles make a non-degenerate isosceles triangle whose three angles are given 
by cos 6*1 = cos 02 — {'ni + l)/4 and cos6'3 = 1 — (to + 1)^/8. The other bounds 
correspond to configurations where the particles are aligned. For to ^ 1, 
the bounds saturates to —9/4 < Eq < —5/4 which is quite rough compared 
to the numerical value Eq > —1.8 obtained for to = 20 with the optimized 
variational method; yet it is better than the results given by the improved 
(Hall-Post) variational method Table 2] with which it coincides for to ^ 1. 
For small to both upper and lower bounds tend to the 2-body exact energy 
and provide acceptable bounds: For instance, when to = 0.05, the differen- 
tial method gives -.34922 s$ £^0 ^ -.29989 while the other ones Table 2] 
give —.59525 ^ Eq (naive variational method) , —.34922 ^ Eq (improved varia- 
tional method), —.34666 ^ Eq (optimized variational method) and Eq ^ —.3375 
(variational with hyperspherical expansion up to L = 8). 

As already mentioned, the differential method, though being less precise for 
N = 3 than the improved or hyperspherical variational approaches, has several 
advantages: it is much simpler, it provides analytic upper and lower bounds that 
furnish an explicit estimation of the errors and, at last but not least, can be 
easily extended to larger N (see below) ; though possible in principle, the gener- 
alization of the improved variational method has not been done beyond A'^ = 4. 



3.3 Several examples of four-body systems 

The optimized variational method has been successfully proposed for A = 4 in 
0] for potentials with scaling-law behavior. For Coulombian interactions with 
a common coupling constant set to —1, tables and figures [31 and |^ compare the 
variational results to those obtained from the differential method when D = 3. 
The same conclusion as in the previous section can be drawn and here are some 
examples of explicit analytic bounds that are obtained by partitioning Fi in 
subclusters made of 3 particles: 

Let us take toi — to,2 — — 1 and toq = to. We have 

9 3m(m^ + 6m + 13) ^ ,„„ , 

O^m^l; \ L^Eo; 32a 

8 8(to +1)^ 
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m 


Naive 


Hall-Post 


Optimized 


Variational 






^(d.„.,_^^(d.„,., 


0.05 


-0.59525 


-0.34922 


-0.34666 


-0.3375 


-0.3246 


0.02467 


-0.3492 


0.1 


-0.6818 


-0.436055 


-0.43434 


-0.423465 


-0.3926 


0.04344 


-0.4361 


0.2 


-0.8333 


-0.58055 


-0.58045 


-0.55915 


-0.5125 


0.06806 


-0.5806 


0.5 


-1.16667 


-0.86805 


-0.86705 


-0.8242 


-0.7813 


0.08681 


-0.8681 


1 


-1.5 


-1.125 


-1.125 


-1.067 


-1.0625 


0.06250 


-1.1250 


2 


-1.83333 


-1.3889 


-1.37135 


-1.30225 


-1.2639 


0.12500 


-1.3889 


5 


-2.16667 


-1.8472 


-1.61705 


-1.53935 


-1.5000 


0.27778 


-1.7778 


10 


-2.3182 


-2.49175 


-1.731 


-1.6495 


-1.6136 


0.37190 


-1.9855 


20 


-2.40475 


-3.74775 


-1.7972 


-1.7134 


-1.6786 


0.43084 


-2.1094 



-6 




0.1 1.0 10.0 m 



Figure 4: Estimations for the ground-state energy Eq of an attractive Coulom- 
bian system made of three particles with masses (m, 1, 1). Naive (A), HaU-Post 
(also called improved, dashed-dotted link) and Optimized (x) provide lower 
bounds on Eq while Variational (open circles with a dashed link) provides an 
upper bound. The corresponding data are taken fromjjj. In the graph, the 
value of iJg^'™ ' has been subtracted. The error-bars are centered at zero (filled- 
circles) and have a width 2AEq'^''^^'\ 



1 < TO < 2V3 - 1; 



9 3TO(TO2 + 6m+13) 
8 8(to4- 1)2 



5m^ + 13to + 2 
2(m+ 1)2 



(32b) 
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2\/3 - 1 to; 

9 3to(to2 + 6to + 13) 3to(to + 1 - 2\/3)2 Sto^ + 13to + 1 

"8 8(m + 1)2 + 8(m+l)2 =^ -^o ^ 2(to + 1)^ ' 

(32c) 

The configuration tliat minimizes i?|^j(g4) given by p5|l corresponds to a 
tetrahedron with an equilateral basis made by particles 1, 2 and 3. The three 
other faces, with particle at one vertex, are identical isosceles triangles, namely 
those which maximize (130(1 when to ^ 3. Such a tetrahedron can indeed be 
constructed provided the angles at particle are lower than 27r/3 which re- 
quires TO ^ 2\/3— 1 ~ 2.464. For to ^ 2\/3 — 1, the configuration that mini- 
mizes £'|^j(g4) seems numerically to correspond to an equilateral triangle made 
by particles 1, 2 and 3 with particle at its center (the flat tetrahedron obtained 
when to = 2\/3 — 1). This simple configuration allows to conjecture the analytic 
lower bound in (|32c|) . For to ^ 1, the upper bound is obtained when the four 
particles are aligned with the 0**^ at one extremity. For to ^ 1, I am not able to 
propose an analytic expression for the upper bound. 

For TO2 — TO3 — 1 and toq = toi = to, a tetrahedron that maximizes all 
the contributions of its faces simultaneously can be constructed for m^^ ^ 

TO < TO^ (-1 + \/T7+ VU- 2VT7)/2. Two identical faces (see figure dl 
corresponding to particles with masses (1,1, to) have their three angles given 
by cos6'i = COS02 — {m + l)/4 and cos^s — 1 — (to -|- 1)^/8. The angles 
{9'i — ^21^3) of the two other faces corresponding to particles with masses 
(to, to, 1) are obtained replacing to by in the previous expressions. For 
1/3 < TO < 3 such faces can indeed be constructed but the pairs of identical faces 
can be put together to construct one tetrahedron provided only that 63 ^ 29[. 
This last condition leads to + 2rn? — lArri? -\-2m-\-l ^ 0. to* and are the 
two positive roots of the four-degree-polynomial, the two others being negative. 
We get 

TO"^ ~ .3622 < TO < TO* ~ 2.7609: ^ — sC . (33) 

2(to + 1 



3.4 Arbitrary number of identical particles plus one dif- 
ferent from the others 

As we have already seen for identical particles, the possibility of partitioning 
the local energy in contributions involving more than two particles allow to 
get some bounds on the ground-state energy of systems made of an arbitrary 
number N of particles. For N ^ 5, this is beyond the scope of the existing 
optimized variational methods. To see one more last example, let us generalize 
both cases of sections 13.11 and 13.21 and consider a system made of one particle 
with mass toq — m and A'^ — 1 identical particles of mass toi = • • • = toat-i = 1. 
All the identical particles interact with the same coupling constants Cij ~ —1 
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-2.29455 


-1.17301 


-1.17283 
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-1.10 


0.07 


-1.1730 


0.1 


-2.65909 


-1.54679 


-1.54167 


-1.45802 
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0.16 


-1.5468 
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-2.47917 


-2.47618 


-2.28857 
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0.28 


-2.4792 


1 


-4.5 


-3 


-3 


-2.78762 


-2.7500 


0.2500 


-3.0000 


3 


-5.625 


-3.9375 


-3.73167 


-3.45553 


-3.3024 


0.6149 


-3.9173 
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-6.3409 


-6.48655 


-4.20877 


-3.90826 


-3.6691 


1.0575 


-4.7266 
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-6.70545 


-40.1395 


-4.4673 


-4.154310 


-3.8412 


1.3266 


-5.1678 


500 


-6.741 


-190.128 


-4.49338 


-4.17914 


-3.8574 


1.3545 


-5.2119 


+00 


-6.75 


— 00 


-4.5 


-4.19259 


-3.8615 


1.3615 


-5.2231 




100.0 m 



Figure 5: Same conventions as in figure 01 for an attractive Coulombian system 
made of four particles with masses (m, 1, 1, 1). 

and the 0'^ interact with eo,i — —.9^- The local energy is given by (|15|l as 



1 



(d - 1)2 I m + 1 



{N-l) + Un - 1){N - 2) + FN-l{qN-l) 



N-l 



a{m ^) ^ |^cos(6?7) -f cos(C)7i) + a(m) cos(iO}) \ (34) 



ij=i 
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0.0204 
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-2.06426 
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0.2986 
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-2.67552 
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0.2744 


-2.6778+ 
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Figure 6: Same conventions as in figure^ for an attractive Coulombian system 
made of four particles with masses (m,m, 1, 1). Numbers signaled with a ^ are 
the lower bounds H33|l . 



where now qn-i stands for the configuration of the A'^ — 1 identical particles. 

a(m) 2g'^ / {m + 1). By simultaneously bounding the contributions of the 
triangles that include the O*'' particle and the contribution of the remaining 
cluster of the identical particles, we immediately have analytic expressions for 
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m 




1 



O3 



m 



Figure 7: For a special range of m, a four-body configuration with masses 
(m, m, 1, 1) that minimizes the local energy can be constructed from the optimal 
configurations of each of its three-body subsystems. 

bounds on Eg. For instance the lower bound is given by 



where F^'^^ is given by and supgj^_^ -Fjv-i has been explicitly estimated 
in section im 

4 Arbitrary two-body interaction 

4.1 Behavior at large distances 

Considering attractive Coulombian interactions is relevant for heavy quarks 
models at short distances but, of course, other kinds of effective potentials are 
required in most models. Since in general no analytic expressions are known 
for the two-body ground-state energies 6^, no method is expected to provide 
explicit non-trivial bounds on Eq. However, if one has some experimental clues 
about eij (by measuring 2-body masses or dissociation energies) or numerical 
estimates as well, it is always interesting to obtain some relations between the e's 





)) (35) 
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and the ground-state energies of larger systems. As mentioned in the introduc- 
tion, this have been achieved in |2] for = 3 and in 1^ for iV = 4 when the 
interactions are of the form Vij{r) cx sign(/3)r'^. 

For /3 > 0, the semiclassical argument given at the end of section |2 shows 
that ()13|l is expected to be unbounded; then 10} gives no information. If we want 
to take the advantage of the simple form H13|) (that is, to keep the choice (O 
with (|ll|l for the test functions), we have to work with finite range potentials. 
When at large distances, the potential is still confining (/3 = 2 for harmonic 
forces 01 [3 — \ for interquark force in quantum chromodynamics P| and [Sj 
for an up-to-date review]), some different ansatz for Lp must be constructed 
in order to eliminate the singular behavior of the ti^'s at infinite distances. 
Actually, the Coulombian case considered in the previous section can be seen as 
an example of a problem where simple poles at finite distances can be eliminated. 
Anyway, in many situations, the 2-body ground-state is expected to depend on 
the behavior of the potential at large distances by exponentially small terms 
only. If, in the integral Q, we decide to keep only those configurations qM 
whose size remains in a physical domain bounded by a cut-off length A, then we 
expect to make an exponentially small error on the estimates of Eq\ this is due 
to the exponential decay of $o when two or more particles separate off. Like 
the Cy's, A is typically obtained from a 2-body dynamics but its precise value 
is irrelevant if the extremal values of the local energy do not depend on it. It 
is precisely the case of the Coulombian interactions (more exactly, interactions 
that can be modeled by Coulombian potential at the energy scale where the 
ground-state exists) for which the local energy (|15|l is invariant under dilations. 

Since, in the present section, we just want to sketch some main guidelines 
without working through the details neither being exhaustive, we will consider 
only the cases where 



4.2 Fitting the 2-body ground-state wavefunction 

What is new, here, is that the differential method allows us to choose directly the 
2-body ground-state wave- functions, or rather their logarithms Sij . Once some 
numerical estimate of eij is obtained in one way or another, we can completely 
bypass the problem of modeling the 2-body potential uniformely. Being free of 
any integration, the differential method can deal with rather complicated, and 
therefore rather realistic two-body test functions. An explicit choice of S'^'s 
provides an explicit form for the local energy (|13|l . 

It frequently happens that we know from experiments the behavior of the 
two-body potential in some specific regimes (most generally at short and large 
distances) but not uniformly. We can therefore, in each of these regimes, tenta- 
tively obtain, with the help of the differential equation Hll|l , the local functional 
form of the two-body ground-state wave- function. Matching these local solu- 
tions together, and then dealing with quite complicated global expression for S 
and w, do not represent a serious obstacle for the computation of (fT!^ . 
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To be a little less speculative, let us consider A'' identical particles with unit 
mass, interacting with a two-body radial potential v{r) such that v{r) — > when 
r ^ oo and 

v{r) — > vor>^ (37) 

r— >0+ 

for some vq and with considering only one case, say /? > 0. The two-body 
stationary Schrodinger equation becomes 



1 d 

dr^ r dr 



2 — ) '^ + ^'^ = £00 (38) 



where eo < will denote an estimate of the two-body ground-state energy; it can 
be considered as another parameter that should fit the experiments involving 

two bodies. Let us guess the behavior of S{r) =^ ln0(r) at short distances by 
writing for 7^ — 1; 

S{r) — > r'^+i . (39) 

^ ' r^0+ (7+1 ^ ^ 

Identifying the leading orders after having reported (|39(l in (|38|l , we necessarily 
get (for D > 1) CT = 1 and sq = — eo/o. The next term in the development of S 
can also be determined. For < /3 < 2, it depends only on the leading term 
(EH) and we have 

This local asymptotic series must be matched with the semiclassical behavior 
at large r 

S{r) -V^r (41) 

r — >-\-oo 

since we have supposed (I36II . The additive constant in S is irrelevant since the 
local energy does not depend on the normalization of 0. A simple choice that 
ensures the local energy to remain uniformly bounded is to take for S' a fraction 
like 

= ■ m 

Figure |H1 show the corresponding for an arbitrary choice of parameters 
together with the corresponding v whose complicated analytic expression on 
[0, -|-oo[ is not needed. 

4.3 Crude bounds 

From equation (|13|l . an immediate upper bound on Eq is given by: 

^N{N - 1) £0 - In{N ~1){N- 2) ^ Eo 

^lN(N~l)eo + ^N{N-l){N-2)a^ (43) 
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-2.546 



Figure 8: D = 3, eo = -1, wq = 1, /3 = 1-5, the function S' given by 
is plotted as well as the corresponding wavefunction = exp S which is the 
ground-state of the potential v for energy eg- The latter can be reconstructed 
from (p with the Schrodinger equation. 

where 

a =^ sup \S'\. (44) 

[0,+oo[ 

Because the constraints between the angles j,i,k are not taken into account, 
these inequalities are expected to be rather rough and their quality deteriorate 
for large N: when positive, the upper bound becomes irrelevant since we already 
know that Eq ^ 0. Indeed, the decreasing of the term with N does not 
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guarantee that the lower bound is better than the Hall-Post bound N{N—l)eo/2 
or even than the naive one N{N — l)eo/2 |3I § 2.1 and 2.2] where eo (resp. eo 
and eo) is the ground-state energy for a particle of mass {N ~ l)/2 (resp. A''/4 
and 1/2) in the central potential v (recall eg < eo < eo)- 

4.4 Reduction to a finite number of Coulombian cases 

In fact, we can find some bounds of (|13|l by reducing the problem to a finite 
number of Coulombian-like cases, that is, where the function to be bound in- 
volves constant factors in front of the cosine (compare H15|l to l|13() ). To see this, 
split the coordinates qi^ into a scaling factor A ^ and some angle variables 9 
among which [N — 1)d — 1 are independent. Each distance writes r^j — \pij{6) 
where the 's are functions that do not depend on the global size of the con- 
figuration but on its shape only. Now, from (|13|) . we define (recall = 1) 

GN{qN)= S'{n,)S'{r,k)cos{jXk) (45) 

and we have 

inf (Gjv(gjv)) =inf (Gjv(0)) (46) 

9n 6 

where 

GN{e) = inf J2 ^'[Vy (^?)] S'iXpMO)] cos(jXfc) • (47) 

Analogous relations are obtained for the maxima. For fixed 9, when A varies 
from to -|~oo, the map A i-^ ('^'[-^Py (^)])o<i<j<Ar_i defines a curve Co in a n- 
dimensional space with n = N{N — l)/2. Cg is bounded if S' is bounded. More 
precisely, Ce is inside the n-dimensional hypercube B '= [(Tmin, Cmax]" where 

a^in= inf S' (48a) 

[0, + oo[ 

and 

CTmax = sup S' . (48b) 

[0, + oo[ 

If S" has the form shown in figure|Sl Ce starts at the origin (S"(0) = 0) and ends 
at the point (—1, . . . ,—1). Taking all the points in B rather than the points 
of Cg leads to a lower bound of Gjy: 

inf Sij SikCOs{j,i,k) Gn{9) (49) 

where s = (s^ )o^i<j^A'-i. Now, whatever the values of the cosines may be, 
the quadratic function in s appearing in the left-hand side of H49|l reaches its 



22 



minimum at a vertex of S^. Let us denote by V, the finite set of the 2" vertices 
s of S {i.e. for all s in V, each Sij is either CTmin or (Tmax)- We have 

inf [inf {fP{0))] < inf (Gjv(gjv)) ■ (50) 
sev e ^ gjv 

where ^ 

/■'v ■ ^' ^ SyS^fcCOs(j,i,A:) . (51) 

In fact, what we have done by obtaining the left-hand side of (|50() is to make the 
values of S"(ry ) independent from those of 9. It follows that the inequality H5U|I 
will be strict if the value of s at the minimizing vertex are incompatible with the 
geometrical constraints on the configuration of the N points. We will illustrate 
this point in the next subsection. We have obtained 

Eo s; l-N(N - l)eo - inf [inf (Ff ^(0))] . (52a) 



and similarly 



In{N - l)eo - sup [sup (Ff ^(f?))] s: Eo . (52b) 



The function Fj^ \0) has the same form as the second sum of the right-hand 
of (|15|1 . Therefore, the computation of the bounds in H52|l is equivalent to a finite 
number of Coulombian problems (with not necessarily attractive interactions 
since the sign of Sij may change) where we must consider all the possible s 
whose components are either (Tmin or CTmax- 



4.5 Three bodies 

As we have seen in section |31 even in the purely attractive Coulombian cases, 
an analytic expression of the extrema of Fn is not known in general. Anyway, 
one can always group in clusters the terms involved in (|51|l like in H24|l . then use 
inequalities like (|25|l and reduce the number of particles. Let us then consider 

^For any constant A, any n-dimensional vector B and any symmetric n X n-matrix C with 
vanisliing diagonal coefficients, the critical points of fn{s) = s- Cs + B ■ s + A are always 
saddle points: the direction S2 = itsign(Ci2)si and Si = for i > 2 makes / increase/decrease 
as ±|Ci2|sj from its critical value. Therefore the extrema of / are reached on the boundary 
of the domain of s. For s restricted to a n-dimensional squared box whose faces are given 
by fixing one Si, the restriction of /„ to one face, i.e, to a (n — l)-dimensional box, leads to 
a function fn—i to which the above argument may be applied again. By repetition down to 
n = 1, we see that the maximum and the minimum of / is necessarily reached at one of the 
vertices of the original n-box. 
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TV = 3. It can be shown® that 



inf [ inf (9))] = 2a^i„a„,ax - ci^ (53) 

with definitions H44|) and H48|l . infg F^^\6) is obtained for an ahgned configura- 
tion which is genericaUy incompatible with s being a vertex of the cube [cTmin, fmax]^- 
For instance, suppose that S' has the shape depicted in figure |H1 where a = 
-CTmin > cTmax > 0; the valuc 2crminO'max-a-^in IS obtained for s = (soi, so2, S12) = 
(cmin, CTmin, CTmax) and should be reahzed for tqi > 1, ro2 > 1 and ri2 ~ rmax 
(the unique finite distance at which S" vanishes); but this is incompatible with 
the alignment condition cos(l, 0, 2) = —1 where particle is in between the two 
others which implies ri2 = rpi + ro2. The inequality (|50|l is therefore strongly 
strict. It can be improved by reducing the size of the cube B to make its mini- 
mizing vertices compatible with the aligned configuration. It can be shown that 
for 5" of the form shown in figure |S1 we have 

2amax5'(2ro) - ct^x «S inf (63(93)) (54) 

where tq is the unique strictly positive distance where S' vanishes^. Even though 
the inequality is still strict because ri2 — 2ro 7^ tqi -I- = 2rmax in general, 
the bound is much better than (|53|l . For instance, if we take the value of the 
parameters corresponding to figure |S| we have 

2CTmi„CTmax " ^1;^ ^ -1.2705 < 2CTmax^'(2ro) - CT^,, ~ -.205 . (55) 

to be compared with the result of the numerical minimization of G3 

inf (63(93)) --.1150 
93 

obtained for the aligned configuration where rpi ro2 
The other bound 

snpiG^iq^)) =snp[snp{FP{0))] 
93 sev 

is actually obtained at the vertex s — (cr, a, a) for an equilateral configuration 
where the common distance rgi = ro2 — 7'i2 is where |iS"| reaches its maximum. 

* The extrema of F^" ^ (6) when 6 varies can be calculated with the help of the appendix with 
a3=soiSo2, 12 = so2Si2 and ai =soiSi2. From ED, we get /(ai,a2,a3) = | (sqi + «02 + *12) 
which is always positive and larger than ai + a2 — ag, —ai + a2 + and ai — a2 + 03 
corresponding to the three aligned configurations for different ordering of the particles. FYom 
I62i . the minimum of Fg''^(9) must therefore correspond to an aligned configuration. Its 
maximum is reached for the configuration described just after equation I57i . 

^Any aligned configuration with roi ~ ro2 ^ '"0 and ri2 = roi -I- ro2 > ''0 corresponds 
to a negative G3. Therefore, as far as its minimum is concerned, the configurations leading 
to a positive infgFg"' can be forgotten (see note 151 . It is straightforward to check that 
all the possible relative positions of Vij with respect to rmax and rg that are compatible 
with ri2 = roi -I- ro2 provide a G3 such that I54i . 



(56) 



ri2/2~ .6107. 



(57) 
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For 5" like in figure|Hl it corresponds to a very large triangle (roi ^1) where a ~ 
1. For D = 3, eo = -1, = 1, /9 = 1.5, from and results lEnj, (|57I), the 
inequalities ifS^ give 

-3 - 3/2 = -4.5 sc £:o -3 + .1150 = -2.885 (58) 

corresponding to a relative error of A£;^'^ '" V-E^o'^ "' '* - 22%. This not really 
impressive (again we emphasize that we are not looking for numerical perfor- 
mance at this stage of development of the differential method) but can be seen 
as an encouraging starting point since the interactions involved so far in the 
three-body system are highly not trivial. It would have required much more 
numerical work to obtain a rigorous window for with variational methods 
(specially a lower bound since the potential considered here does not follow a 
power-law behavior). 

5 Conclusion 

The differential method appears to offer a completely new strategy for estimat- 
ing a ground-state energy. For many-body systems, we have seen on several ex- 
amples how this approach can be fruitful. For attractive Coulombian particles, 
it can compete with existing others methods (that are based on the variational 
principle) on several levels: it provides upper and lower bounds with compara- 
ble numerical precision, its simplicity renders the analytic calculations tractable 
even for large N and/or allows a low cost of numerical computation. Beyond 
purely Coulombian systems, the differential method, being so general, offers a 
remarkable flexibility. As have been sketched in the previous section, one can 
deal with systems where interactions can be very rich (possibly short-ranged 
with an a priori cut-off); several regimes which are valid at different scales can 
be implemented at once. There is some hope that future works successfully 
apply the differential methods for proper realistic potentials. 

Unfortunately, I have not been able to generalize the differential method to 
fermionic systems where the ground-state spatial wave-function is antisymmet- 
ric. In such cases, the presence of non-trivial nodal lines (Sj breaks down the 
proof of inequalities Q . 

There is a lot of work to be done regarding a systematic improvement of the 
bounds, once some finite ones have been found with a given at first attempt. 
In this paper, we have not considered some free parameters on which a (fam- 
ily) of test functions, say (^<^, may depend. As shown for a one-dimensional 
system |15| . the locality of the differential method may require a very few 
number of A's at each optimization step (unlike for the variational test func- 
tions) for obtaining substantial improvements of the bounds by calculating, say 

supj [infq (-t'loc ■ ^ precise proof that this approach is efficient for several 
dimensions remains an open interesting problem. 

I thank Jean-Marc Richard for a critical reading of the first proof of this 
manuscript and acknowledge the generous hospitality of Dominique Delande 
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and Benoit Gremaud of the group "Dynamique des systemes coulombiens" at 
the Laboratoh'e Kastler-Brossel. 



Appendix: Extrema for the three-body Coulom- 
bian problem 

For the three-body Coulombian problem, as can be seen from the second sum 
in the right-hand side of 1)15(1 . we must find 

F™''''(ai, 02, 03) =^ sup [oi cos 01 -I- 02 cos 6*2 -I- 03 cos ^3] (59) 

triangles 



and 



def 



3 (01,02,03) = 



inf [oi cos 6*1 -I- 02 cos 02 + 03 cos 6*3 



triangles 



(60) 



where the ^'s are the angles at tlie three vertices of the triangle made of the 
three particles. The o's are some real parameters that depend on the masses 
and the coupling constants. Let us define 



/(oi,a2,a3) = - 



ai02 , 0103 



02 



01 



(61) 



then we have 

max 

F^™ (oi, 02, 03) = ™t^{ai-|-a2-03, 01-02-I-03, -01-1-02-1-03, /(oi, 02, 03)} (62) 

/ is considered in the list (|62|l if and only if the following three conditions are 
satisfied simultaneously 



02 
03 



as 

02 



0203 



s$ 1 



Ol 
02 



02 
Ol 



O1O2 



OS 



^ 1 ; 



1 



03 

Ol 



Ol 

03 



0103 



s$ 1 



(63) 



Here is the proof: We will restrict the values of the 0's to [0, tt] and the con- 
straint 01-1-02 + ^3 = TT is implemented by the Lagrange multiplier method. 
We are led to extremalise the function 03(61, 62, 83) = aicos0i -I- 02 cos 02 -I- 
a3cos03 -I- £(01 + 02 + 03 — tt) for unconstrained (0i,02,03) € [0,7r]'^, £ being 
the Lagrange multiplier. The three conditions do^G3 = for i = 1,2,3 lead 
io i = oisin0i = O2sin02 = O3sin03. The case £ = corresponds to the 
alignment of the three particles and gives the three first values in the list H62|l 
corresponding to (0i,02,03) = (0,0, tt) and its circular permutations. 

Taking into account the constraint on the angles, we have H — O3sin03 = 
03 sin 01 cos 02-I-O3 sin 02 cos 0i — £(03 cos 02/01-1-03 cos 01/O2). Therefore when i ^ 
0, we find oicos0i -I- O2cos02 — 01O2/03 as well as the other relations that 
are obtained by circular permutations of the indices. From the decomposition 

F3 = i(oi cos 01 -I- O2 cos 02) -I- i(02 cos 02 -|- O3 COS 03) + i(oi COS 0i -f O3 COS 03), 
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we obtain the value that must be considered in (|62|l if and only if there 
exists some 0's such that 

fli sin 01 = a2 sin02 = ^3 sin^a and di + 62 + O3 = tt . (64) 

Solving these three equations leads to the values for the three | cos 9i \ that are 
precisely given by the left-hand sides of the inequalities Ht)3|l . 
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